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ABSTRACT 



Some experiments on LEAR obtained unusual behavior of the pp interaction near 
the threshold. The experiments on pp forvard scattering detected zeros and big 
variation of p and at the same time a smooth rising of a to t with lowing energy. 
Many models has difficulties in explanating this fact. In the PS170 experiment with 
a good statistical accuracy the unexpected behavior of the proton electromagnetic 
form factor was found. All these experiments can be considered as an indication for 
the existence of a low lying pp bound state 'baryonium'. This statement coincides 
with that made for interpretation of the energy dependence of the total cross-section 
e + e~ — > hadrons in FENICE. There is a model (based on analyticity) which explained 
aforementioned experiments and the fact that this 'baryonium' is not seen in the 
OBELIX pp annihilation cross-section. Thus LEAR experiments and FENICE one 
are consistent near pp threshold and compatible with the existence of 'baryonium'. 

PACS:13.40.Fn — Electomagnetic form factor; electric and magnetic moments. 
14.20 — Baryons and baryons resonans(including antiparticles). 
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1 The database and previous knowledge 

The experiment on LEAR which is a part of the CERN antiproton complex gives 
a rich information about low energy antiproton physics. The experiments (PS172, 
PS173) HI], H on pp scattering give the data on da/dQ, a tot and p. To search for 
bound state cross section measurements are the most straightforward experiments 
to perform. The analysis of da/dQ gives an indication of the bound states near pp 
threshold 0. Some of them are consistent with the strong interaction shifts and 
width of protonium |3J]. A resonance (a bound state having a mass bigger the pp 
threshold) may be seen as a bump in a tot . But the measurements of the pp total 
cross section above 180 MeV/c indicate its smoothly varying behaviour 0. The 
most remarkable result in pp elastic scattering has appeared in the data on real-to- 
imaginary ratio of the forward scattering amplitude p which at LEAR was measured 
down to 180 MeV/c ||. For the range 350 <pi < 700 MeV/C the behaviour of p can 
be explained by insertion of a pole below threshold in the dispersion relation analysis 
0. But LEAR measurements 0, below 350 MeV/c indicate that the p is turning 
upward again. The reason for this unusual behaviour is not yet clear. It might be 
caused by a pp bound state M but not by an nn threshold M. Experimental the p 
was always determined from elastic differential cross section in the Coulomb-nuclear 
interference region. The method used to extract p from such data sometimes has 
been criticized ||. But at high energies the method is consistent with the predictions 
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of dispersion relations. So p from [Q, |J will be considered below as reliable. 
The results of experiment PS-170 on the study of annihilation pp — s> ee at low 
energies [10] have no adequate interpretation till the present day. They resulted 
in an unexpected behaviour of the proton electromagnetic form factor near the pp- 
threshold in the time-like region, where s < 4.2 GeV 2 . The data on 
| G \—\ G m>p | = | G £iP | point to a large negative derivative at the threshold that 
rapidly grows to zero or even to positive values at s ~ 4 GeV 2 . The magnitude of 
the derivative at the threshold is determined by the threshold value | G \ — 0.53 ± 
0.05. One of the early values, | G |= 0.51 ± 0.08, does not contradict the results 
of ref. [10]. It was obtained [11] from the ratio of frequencies of pp annihilations 
at rest into ee and 7r + 7r~ pairs in liquid hydrogen. The determination of | G | at 
the threshold is a complicated problem since one should simultaneously consider the 
Coulomb and strong interactions in the pp-system and requires some approximations. 
These approximations have been analysed in ref. [12] where a new scheme is proposed 
for the determination of | G |. This scheme gave the value | G \= 1.1 that confirms 
the results of ref. [10]. Quite recently, a new attempt has been undertaken for 
determining | G | at the threshold [13]. Combining the data on widths of pp-atoms 
obtained in the synchrotron trap with the results on the low-energy annihilation cross 
section in pp-system, the authors concluded that | G \= 0.39 or even | G | = 0.3. This 
allows us to infer that there is no abrupt change of | G | at the threshold. Thus, 
the authors of [12,13] propose a new view on the method of calculating | G | at the 
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threshold from experimental data. 

Let us now proceed to works that suggest the interpretation of the results of the 
experiment [10]. In ref. [14] an attempt is made to consider the interaction in the 
final state. The basic result is the formula G = ce lS , where c is a slow variable 
function of q 2 at the threshold (q is the momentum in c.m.s. of the pp-system) and S 
is the NN scattering phase. Since the phase 5 is complex at the threshold, we have 



where a is the complex scattering length. Owing to | G | being linear in q, the 
quantity d \ G \ /ds is infinite at the threshold. Analysis of the first four points from 
[10] with respect to the x 2 -criterion gives the values: | c |= 0.53±0.02, Im a = 0.62± 
0.08 fm, x 2 = 0.07. The authors of [14] employ the values: | c |= 0.52, Im a = 0.8 fm; 
they identify Ima with the quantity Im a( 3 5i) computed from the experiment [15]. 
The description is qualitative since x 2 ~ 10- The authors of [16] assert that a good 
description of all known data on nucleon electromagnetic form factors, including the 
data of [10], is obtained on the basis of a new formulation of the vector-dominance 
model (VDM) and its subsequent unitarization. In what follows, we will use different 
models of that type, therefore we consider them in detail. They are based on the 
expressions for the Dirac and Pauli nucleon form factors in VDM: 



G 



—\ c 



1 — q • Im a 



(1) 



Fn{s) = E 



V 



fv,NN m 2 
f v ml- s ' 



(2) 
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where m v is the mass of a vector meson, f V NN is the coupling constant of a vector 
meson with a nucleon, /„ is the universal constant in the so-called identity of current 
and field. Imposing constraints on the parameters of formula (2), one can easily 
find the experimental value Fn(s = 0) and the asymptotics following from the quark 
counting rules [17] that coincides with the QCD-asymptotics within the logarithmic 
accuracy. Then, the model is unitarized with the help of a uniformizing variable. 
As a result, vector mesons acquire widths, and the form factors can be calculated 
for all s. So, all experimental data can be described both in the space-like (s < 0) 
and time-like (s > 0) regions. Satisfactory description of more than three hundred 
values of | F N | requires about ten free parameters in the formula (2). Besides, 
this approach allows a model-dependent reproduction of the form of ImF N , ReF N 
in the whole time-like region. This fact will be used below. Results of the analysis 
according to this scheme are presented in ref. [16]. The data of the experiment 
PS-170 are explained by including the third radial excitation p(770) with the mass 
y/s = 2.15 GeV into formula (2) and are plotted in Fig.l. 

2 Formulation of the analytical model 

It is easy to see that the nucleon form factor, according to formula (2), has the 
following imaginary part 



lmF N = iiY, 



V 



m: 




5(s-m 2 v ). 



(3) 
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Formula (3) is an approximate expression obtained from the unitarity condition which 
allows one to reproduce equation (2) with the use of dispersion relations for F N . We 
write the starting expression for the unitarity condition as follows: 

Im < o | | NN >=J2< ° \ J^\n><n \ T + \NN >, (4) 

n 

where is the electromagnetic current of a nucleon N, and \n > is the complete set 
of admissible intermediate states. In our case, it is of the form 

\n >= \2n >, \3tt >,..., \KK >, \N, N > . (5) 

Frazer and Fulco [18] were the first who computed the contribution of the two-pion 
state and predicted the p-meson on the basis of data on F N . By choosing different 
terms in the sequence (5), one can obtain many models of the type (2). Earlier, the 
model of ref . [19] was used in [20] and the contribution of an NN intermediate state 
was calculated. This contribution is important for two reasons. First its consideration 
results in a new branch point in formula (2), the threshold of the reaction NN situated 
on the lower edge of the energy region studied in ref. [10]. Second bound states or 
resonances in NN -system near the threshold will influence the behaviour of F N (s) 
in the nonobservable region below the iViV-threshold and in the observable region 
above the iViV-thresfiold investigated in ref. [10]. It is clear that the state |iViV > 
appears on the background of the sum of other states of the series (5) and the result is 
model-dependent. Therefore, it is important to study the degree of that dependence 
by considering another model differing from the one used in [20] for F N (s) as a 
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background for the state \NN >. 

We will take the model of ref. [21] formulated in terms of the Sachs form factors G 
measured experimentally. The model is based on the formulae 

G m , P (s)=i: — 6kis \ — , G e , P (s)=i: — (3k(s) / — , (6) 

where e k (s) = 4 + ^ , (5 k {s) = Pl + ®* ■ (7) 

s-a k - 7fcVSfc -s s-ctk- 1kV s k - s 

The energy behaviour of electromagnetic form factors is explained with the use of 

three resonances:p, u, <p specified by indices k — 1,2,3 in formula (6). The masses, 

widths and thresholds a kl ^ kl s k are taken from experiment. The model parameters 

are the coupling constants 

((3{ + e?s)/i(s) = g 1P (s)g pNN (s), 
(I3\ + e° 2 s)f 2 (s) = gju,(s)gu,NN(s), 
(Pi + 4 s )h(s) = ^(s)^jvjv(s) , 

where 

f k (s) = 1 . ■ (8) 

s-a k - 7fcV s fe - s 

This unusual form of the constants is chosen by analogy with the index of refraction in 
optics. They are not only energy-dependent, but also contain a complex component 
when s > s k . The coupling constants are chosen so as to be consistent with the 
known experimental data at s = 0. Then, we are left with two free parameters e° and 
e° to be defined from the conditions required at s — > 00. The SU (3) symmetry should 



hold on the asymptotics identically. This condition seems to be the weakest one since 
it can be changed by including new vector mesons into consideration. Therefore, 
the parameters e° and e° are determined according to the x 2 criterion on the basis of 
experimental points \G P \ cited in refs. [22]. An interesting feature of the model [21] is 
that it correctly describes the ratio |Gr p |/|Gr n | above the pp-threshold. More exactly, 
it reproduces the experimental value \G n (s = 4 GeV 2 )\ = 0.42 ± 0.06 (see [23]). The 
model result for \G P \ is drawn in Fig. 2. and e° = — 3.41, e° = 3.23, x 2 — 10.1. 
The influence of the |iViV > contribution to the unitarity condition (4) on \G\ is 
computed in the same way as in refs. [20, 24]. We construct the analytic model for 
the forward elastic scattering amplitude T in terms of the uniformizing variable 



where s is the conventional Mandelstam variable equal to the square of the total 
energy of a pp-system in the c.m.s. in units M p . The variable z contains branch 
points at s = 0; 4 corresponding to the reaction threshold of elastic pp and pp- 
scattering and an effective branch point at s = a corresponding to the nonobservable 
region for the elastic pp- scattering. The threshold of process pp — > pp is mapped 
into points z = ±1 on the z-plane; whereas the infinit s-plane point, into points 



±zi, ±l/zi, where z\ = y f^f- Disposition of all the four sheets of the Riemann 
surface of the function z(s) is drawn in Fig. 3 for a = 1.44. In ref. [24] it is shown 
that the experimental data on p = ReT/ImT and a to t can be well described provided 



z = 



4(s - a) a(s-A) 



(9) 
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that the pp-system possesses a quasinuclear bound state with the binding energy 
E = (1.88 ± 0.05) MeV and width T = (1.6 ± 0.1) MeV. The scattering amplitude 
was taken in the form 

T Tb+ z-(z p ) 1 z-(z p )2 ^ 

where T b (s) is a polynomial in z, (2,9)1,2 = 1 =F 7 ± i5 and a = 1.44, 10 2 7 = 
—0.54 ± 0.02, 10 2 5 = 2.6 ± 0.08. The pole terms represent the contribution of the 
quasinuclear state; whereas the polynomial determines the contribution of a nonres- 
onance background of S,P and D-waves. Special attention was paid to the threshold 
value of the T amplitude which is complex [24]. The amplitude (10) well describes 
the experimental data up to 4.4 GeV 2 in terms of the variable s. It is valid in the 
vicinity of z — 1 and has two poles in distinction to the usual quantum mechanical 
amplitude. Appearence of the two poles in the variable z inctead of the one pole 
in the variable q in the scattering amplitude T is a consequence of choosing z as 
uniformazing variable. Another important feature of the formula (10) is the form of 
the pole term contribution to Im T and Re T. The bound state (pole) contribution 
to the lmT(pi = 80 MeV/c) is about 10% of the total value Im T pP . On the other 
hand, the bound state contribution to the ReT is larger then the background one 
and ensure the correct value p (see Fig. 4,5). Near the pp-threshold the pole con- 
tribution to the unitarity condition (4) becomes dominant, and thus, we will restrict 
ourselves to the pole approximation. Quantum numbers of this state are unknown. 
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A detailed scheme of calculation corresponds to the scheme by Frazer and Fulco [18] 
for the contribution of different partial waves to lmF N . In our case, it gives that 
these states are either 3 Si or 3 Di. Then, the unitarity condition (4) is reduced to the 
Riemann boundary-value problem [25] that can be solved (see Appendix). Inside the 
ring containing the unit circle (Fig. 3) the solution is of the form 

r n-n 

Gpo1 nU* -(*„)<)(* + (*;)«)' [ } 

where c(z) is an entire function within which the solution is determined. Setting 
c(z) = ci(z) ■ (z 2 — z\)(l — z 2 z 2 )/(l — zf) 2 , we can ensure the asymptotic behaviour 
of Gp i at infinity. Taking advantage of ci(z) being arbitrary, we assume the solution 
to be of the form 

Gpoi(z) 



[z 2 - z\){z 2 z[-\) 



+ 



A \, 1 1 \ _ / 1 1 

+A2 {(^)T + 7=W) ~ (^+tyh + z + {z + {z;))\ (12) 

Around the pp-threshold the equalities | G e . p | = | G mtP \—\G\ hold valid and, under 
this assumption, the experiment in [10] was analysed. Therefore, we put 

G e ,p + G mp = 2G W , (13) 

where the functions Gr e ( m ) jP are given by formulae (6). Considering the contribu- 
tion of the \NN >-state to the unitarity condition (4), we obtain for the proton 
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electromagnetic form factor G: 

G = G w + G pol . (14) 

We shall assume the position of poles to be known from ref. [24]; then, the form 
factor G depends on two free parameters Ai, A 2 . The behaviour of G po i on the upper 
edge of the cut [a, oo) around the iViV-threshold is determined by the poles (z p )i 
and (z p ) 2 ; whereas on the lower edge, by the poles {z*)i and (z*) 2 . If we calculate the 
common denominator for the contributions of the poles (z p )i and (z p ) 2 in the formula 
(12), the energy factor (z — 1) will arise in front of the parameter A 2 ; whereas a 
constant, in front of the parameter A x . This allows us to draw analogy between 
the parameter A 1 and e\,[3\ as well as between A 2 and e° k) f3® in formula (7). The 
expression for G po i (11) follows from the unitarity condition and analytic properties 
of the proton form factor and iViV-scattering amplitude. Therefore, formulae (6) are 
substantiated, irrespective of the above mentioned analogy with optics. The result 
of the analysis (Fig. 4) according eq. (12) is presented in the table 1 and parameters 
are equal: a = 0.23 ±0.04, e° = 2.97±0.03, e° = 3.23, 10 2 A 1 = 0, 10 2 A 2 = 1.2 ±0.01. 

3 Discussion of the results 

The parameters A\ and A 2 representing the coupling constants of a quasinuclear 
bound state are sensitive to the background shape in formula (14) as follows from 
comparision of this fit and the fit of ref. [20] {A x ^ in ref. [20]). The magnitude of the 
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background is determined by the parameters e° and ejj and is slow changing function 
in the s interval under investigation. The parameters A 1 ,A 2 ,a determine the rapid 
change of G in formula (14). Via separating the parameters into these two groups, 
we can obtain their statistically reasonable values (table 1). The analysis would be 
considerably simplified if the experimental values of s > 4M^ were known for ImG 
and ReG. Their determination requires polarization experiments whose theoretical 
study is carried out in ref. [26]. 

Recently two independent experiments gave new information on the pp interaction 
at low energy. The value of the pp annihilation total cross section down to the 
momenta 43 MeV/c have been measured by OBELIX experiment [27] at LEAR and 
no resonant behaviour of the cross section was found. The existence of some stracture 
in the ee — ► hadrons cross section near the PP threshold was indicated in FENICE 
at ADONE [28]. A combined analysis of these data and the data on the proton 
form factor provides a good candidate for the quasinucler bound state with the mass 
M = 1.85 ± 0.01 GeV 2 and the width T = 40 ± 10 MeV. This candidate dosn't 
contradict our candidate [20]. Then the question arise why this candidat is not seen 
in the OBELIX experiment on the pp annihilation cross section at very low energy. 
The first reason for that is the mass of 'baryonium' which is less then 2M p . The second 
is based on our analytical model. In this model (a to t)poie — 0.1 a tot at low energy 
(see Fig. 4) but a ann < a to t- From these inequalities it is clear why 'baryonium' is 
not seen in OBELIX data. On the other hand, in FENICE experiment cross section 
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e + e~ — > hadrons depends not only on ImT but also on ReT for which the pole 
contribution is large. That is the reason why 'baryonium' is not seen in OBELIX and 
seen in FENICE. Thus the results of bouth this experiments are consistent. 
Finaly we mention a pure theoretical result; the method of derivation of formula (11) 
for describing a quasinuclear state can be applied to any vector meson in formula (2). 
Therefore any vector meson will be characterized not only by the mass and width but 
also by two parameters like coupling constants. In other words, the effective coupling 
constants will be energy-dependent, what is assumed in ref. [21] and is reflected in 
formulae (7). 

Appendix 

The unitarity condition (4) is an exact equation if use is made of the complete sys- 
tem of admissible intermediate states (5), otherwise it is an approximate equation 
dependent on the assumptions made. Let us take it in the form 

ImF = F(e iS sin 5)* +g, 

where 5 is the iViV-scattering phase with quantum numbers of the pole state unknown 
yet; g is the contribution of all other processes in the same channel. We reduce it to 
the form 

F = e 2iS F* + 2ig. (A.l) 
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The relation (A.l) is valid for Ims = and Res > AM 2 . The function F is analytic 
in the complex plane s with the cut [4M^, oo) outside of which F*(s) = F(s*). This 
relation represents a linear inhomogeneous Riemann boundary-value problem for the 
function F. If e 2lS has a pole near the cut, then in its vicinity we can consider the 
homogeneous problem 

F = e 2i5 F*. 

As it is known [25] , the main difficulty in solving it consists in constructing a function 
analytic in the plane s and coincident on the cut with e 2%5 . However, if e 2tS is taken in 
the form admitting the analytic continuation onto complex s, the problem is reduced 
to the solution of a functional equation for F in the uniformizing variable z. We will 
represent e 2tS in the form 

>i8_Tr( z - z j)( z + z j) 



=n 



. (z- Zj )(z + z*) 

The function e 2tS is real on the imaginary axis z, i.e. on the real axis s when s < a. 

Equation (A.l) is valid on the cut [4Mp,oo) that transforms into the real axis z = 

x + iy , and F(s) -> F(x), F*{s) -> F(-x) 

(x-z*)(z + Zj ) 
Fix) = — rFl—X), 

where we took only one pole, without loss of generality. The latter functional equation 
for F(x) can be written as follows 

F(x)(x - Zj )(x + z*) = G(x), 
15 



G(x) = G(-x) 
and thus F(z) is representable in the form 

G(z) 



F(z) 



Ujiz-zj^z + z*)' 

where G(z) is an entire even function of the variable z. The inhomogeneous boundary- 
value problem (A.l) can be solved in a similar manner and formula (11) can be proved. 
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Table 1: 



S GeV 2 


G eX p 


\G\ 


Xi 


3.523 


0.53 ±0.05 


0.63 


3.9 


3.553 


0.39 ±0.05 


0.35 


0.63 


3.57 


0.34 ±0.04 


0.32 


0.26 


3.59 


0.31 ±0.03 


0.3 


0.15 


3.76 


0.26 ±0.014 


0.27 


0.66 


3.83 


0.25 ±0.01 


0.27 


1.9 


3.94 


0.247 ±0.014 


0.254 


0.23 


4.18 


0.252 ±0.011 


0.221 


8.1 
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Figure Captions. 

Fig. 1. The curve from Fig. 3 a of ref. [16] on a larger scale. The quality of the fit 

PS-170 data is very poor. 

Fig. 2. Our fit to the old data [22] by G w . 

Fig. 3. Disposition of four sheets of the Riemann surface of the function z(s) 
for a = 1.44. The threshold pp is mapped into points z — ±1. 
Fig. 4. The pole contribution to the ImT. 
Fig. 5. The pole contribution to the p. 

Fig. 6. Our fit to PS-170 data with account of the pole contribution (eq.(12)). 



23 



0,71 




3.5 3.68 3.86 4.04 4.S2 



S GeV2 

Figure 1: 



24 



0,58 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



Q.45B 



O 0.335 



0,313 



Q.09 




3.57 



3.95 



4.33 

S GeV2 



4.72 



5,1 



Figure 2: 



25 




26 



4.2 I 1 1 1 1 1 1 1 1 1 1 r 




200 400 600 



Pi(MeV/c) 



Figure 4: 



27 



0.9 
0-0 



-1.5 - 




-3.9 - 



-6.3 







~i 1 1 1 1 1 1 1 r 



~i r 




□ p, x 2 /ndf 1 

O ppol 



J I I I I I I I L 



200 400 

Pi(Mev/c) 



600 



Figure 5: 



28 




Figure 6: 



29 



